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Abstract 

Using a 2D model of capacitively coupled Josephson junction arrays (created by a network of 
twin boundary dislocations with strain fields acting as an insulating barrier between hole-rich 
domains in underdoped crystals), we study the influence of chemical pressure (V/x) on nonlinear 
(i.e., VT - dependent) thermal conductivity (NLTC) of an intrinsically granular superconductor. 
Quite a substantial enhancement of NLTC is predicted when intrinsic chemoelectric field E^ oc V// 
closely matches the externally produced thermoelectric field Ej- oc VT. The estimates of the model 
parameters suggest a realistic possibility to experimentally monitor this effect in non-stoichiometric 
high-T<7 superconductors. 

PACS numbers: 61.72.-y, 74.25.Fy, 74.50.+r, 74.81.Bd, 74.81.Fa 
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High resolution imaging of the granular structure in underdoped Bi 2 Sr2CaCu20g + s crys- 
tals |l| has recently revealed an apparent segregation of its electronic structure into super- 
conducting domains (of the order of a few nanometers) located in an electronically distinct 
background. In particular, it was found that at low levels of hole doping (5 < 0.2), the 
holes become concentrated at certain hole-rich domains. (In this regard, it is interesting 
to mention a somewhat similar phenomenon of "chemical localization" that takes place in 
materials, composed of atoms of only metallic elements, exhibiting metal-insulator transi- 
tions 2].) Tunneling between such domains leads to intrinsic granular superconductivity 
(GS) in high-T c superconductors (HTS). Probably one of the first examples of GS was ob- 
served in Y BaoCusOjs single crystals in the form of the so-called "fishtail" anomaly of 
magnetization [3J. The granular behavior has been related to the 2D clusters of oxygen 
defects forming twin boundaries (TBs) or dislocation walls within CuO plane that restrict 
supercurrent flow and allow excess flux to enter the crystal. Indeed, there are serious argu- 
ments to consider the TB in HTS as insulating regions of the Josephson SIS-type structure. 
Besides, a destruction of bulk superconductivity in these non-stoichiometric materials with 
increasing the oxygen deficiency parameter 5 was found to follow a classical percolation 
theory 4|. In addition to their importance for understanding the underlying microscopic 
mechanisms governing HTS materials, the above experiments can provide rather versatile 
tools for designing chemically-controlled atomic scale Josephson junctions (JJs) and their 



arrays 
devices 



JAs) with pre-selected properties needed for manufacturing the modern quantum 
5, |f|. To understand how GS manifests itself in non-stoichiometric crystals, let 
us invoke an analogy with the previously discussed dislocation models of grain-boundary 
Josephson junctions (GBJJs) (see, e.g., [7|, |8j and further references therein). Recall that 
under plastic deformation, grain boundaries (GBs) (which are the natural sources of weak 
links in HTS), move rather rapidly via the movement of the grain boundary dislocations 



y, 0, 



ii 



in HTS sin- 



(GBDs) comprising these GBs. At the same time, observed 
gle crystals regular 2D dislocation networks of oxygen depleted regions (generated by the 
dissociation of < 110 > twinning dislocations) with the size do of a few Burgers vectors, 
forming a triangular lattice with a spacing d > d ranging from lOnm to lOOnm, can provide 
quite a realistic possibility for existence of 2D Josephson network within CuO plane. Recall 
furthermore that in a d-w&ve orthorhombic YBCO crystal TBs are represented by tetrag- 
onal regions (in which all dislocations are equally spaced by do and have the same Burgers 
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vector a parallel to y-axis within CuO plane) which produce screened strain fields [lOj 
e(x) = e(0)e~' x !/ d ° with | x | = \Jx 2 + y 2 . Though in Y Ba^CujPi-^ the ordinary oxygen 
diffusion D = D§e~ Vd l kBT is extremely slow even near T c (due to a rather high value of the 
activation energy Ud in these materials, typically Ud — leV), in underdoped crystals (with 
oxygen-induced dislocations) there is a real possibility to facilitate oxygen transport via the 
so-called osmotic (pumping) mechanism jlj . which relates a local value of the chemical 
potential (chemical pressure) M x ) = /x(0) + V/i-x with a local concentration of point defects 
as follows c(x) = e _M ( x )/ fcsT . Indeed, when in such a crystal there exists a nonequilibrium 
concentration of vacancies, dislocation is moved for atomic distance a by adding excess va- 
cancies to the extraplane edge. The produced work is simply equal to the chemical potential 
of added vacancies. What is important, this mechanism allows us to explicitly incorporate 
the oxygen deficiency parameter 8 into our model by relating it to the excess oxygen concen- 
tration of vacancies c v = c(0) as follows 8 = 1 — c v . As a result, the chemical potential of the 
single vacancy reads n v = MO) — — ^B7Tog(l — 8) ~ ksT8. Remarkably, the same osmotic 
mechanism was used by Gurevich and Pashitskii [10] to discuss the modification of oxygen 
vacancies concentration in the presence of the TB strain field. In particular, they argue that 
the change of e(x) under an applied or chemically induced pressure results in a significant 
oxygen redistribution producing a highly inhomogeneous filamentary structure of oxygen- 



deficient nonsuperconducting regions along GB [ll| (for underdoped superconductors, the 
vacancies tend to concentrate in the regions of compressed material). Hence, assuming the 
following connection between the variation of mechanical and chemical properties of planar 
defects, namely M x ) = -^^o^(x) (where Qq is an effective atomic volume of the vacancy 
and K is the bulk elastic modulus), we can study the properties of TB induced JJs under 
intrinsic chemical pressure V/i (created by the variation of the oxygen doping parameter 8). 
More specifically, a single SIS type junction (comprising a Josephson network) is formed 
around TB due to a local depression of the superconducting order parameter A(x) oc e(x) 
over distance do producing thus a weak link with (oxygen deficiency 8 dependent) Joseph- 
son coupling J (5) = e(x)J = J (8)e-W d ° where J (8) = e(0)J = (fi v /Ktt )J . Thus, 
the considered here model indeed describes chemically induced GS in underdoped systems 
(with 8^0) because, in accordance with the observations, for stoichiometric situation 
(when 8 ~ 0), the Josephson coupling J(8) ~ and the system loses its explicitly granular 
signature. 
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There are several approaches for studying the thermal response of JJs and JJAs based 
on phe_lo^of the Josephs™ effect in the presence of thenna. gradients (see, 



e.g 



14, H5|, 
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171 | and further references therein). In this paper, within a 2D model 
of capacitive JJAs (created by a regular 2D network of TB dislocations), we study the influ- 
ence of chemical pressure effects (described by the gradient of the chemical potential V/i) on 
nonlinear (i.e. VT-dependent) thermal conductivity (NLTC) k(T, V/i; VT) of intrinsically 
granular superconductors. As we shall see, in a sharp contrast with its linear counterpart, 
NLTC exhibits quite a pronounced enhancement when the intrinsically induced chemoelec- 
tric field = ^|V/i| matches the externally produced thermoelectric field Et = So|VT|. 

To adequately describe transport properties of the chemically induced granular supercon- 
ductor under a simultaneous influence of intrinsic chemical pressure V/i(x) = KQoVe(x) 
and applied thermal gradient VT, we employ a model of 2D overdamped Josephson junction 



array which is based on the following total Hamiltonian [12] 

H(t)=H T (t)+H c (t), (1) 

where 

N 

n T (t) = J2J t A 1 - cos M t )} ( 2 ) 

ij 

is the well-known tunneling Hamiltonian, and 

H C W = E^ (3) 

accounts for a Coulomb contribution due to mutual capacitances CV, between grains with a 
net charge 

= £ / dt'htf) (4) 

induced by the normal current hjif) = Vij(t)/Rij within the array, where Vij(t) = (f^)^|f- 
and Rij are the corresponding voltage and resistance between grains in their normal state. 

According to the above-mentioned scenario, the tunneling Hamiltonian Tirit) introduces 
a short-range (nearest-neighbor) interaction between iV junctions (which are formed around 
oxygen-rich superconducting areas with phases <pi(t)), arranged in a two-dimensional (2D) 
lattice with coordinates Xi = (xi,yi). The areas are separated by oxygen-poor insulating 
boundaries (created by TB strain fields e(xy)) producing a short-range Josephson coupling 
Jij = J (5)e-^/ d . Besides, as we have seen, the strain fields e(xjj) also control the dielectric 



properties of the array. Namely, the mutual capacitances between grains are defined as CV 



27re e(x ij )|x ij | 



Co(\xij\/d)e l Xi jl/ d . Thus, both the Josephson energy and the charging 



energy of the array vary exponentially with the distance Xy = Xj — x,- between neighboring 



junctions (with d being an average grain size). r . 
coupling is governed by the standard expression 

A(T) 



he temperature dependence of Josephson 



Jy(T) = Jy(0) 



A(0) 



tanh 



A(T) 
2k B T 



(5) 



where Jy(0) = [A(0)/2](R /R i j) with A(T) being the temperature dependent gap param- 
eter; Rq = h/Ae 2 and Rij oc Jy = R n e [yLi3 ^ d stand for quantum and normal resistance, 
respectively. For the explicit temperature dependence of the gap parameter we used the 
analytical approximation 
temperatures. 

By analogy with a constant electric field E, a thermal gradient VT applied to a chemically 



12], namely A(T) = A(0) tanh (2.2 ) which is valid for all 



induced J, 



follows iia, y, lis | 



A will cause a time evolution of the phase difference across insulating barriers as 



2e(E M - E r )x 



-t (6) 
T^V/i and E T = SqVT are the 



Here is the initial phase difference (see below), E^ 
induced chemoelectric and thermoelectric fields, respectively. So is the Seebeck coefficient. 
In what follows, we assume, for simplicity, that V/i = (V^/i, 0, 0) with V x fi = Afi/d, and 
that VT = (V x T,0,0). 

In order to study the most interesting situation when the externally produced thermo- 
electric field Et becomes comparable with intrinsically induced chemoelectric field E M , let 
us consider a nonlinear generalization of the conventional Fourier law and the resulting 
nonlinear thermal conductivity (NLTC) under the influence of chemical pressure. In what 
follows, by the NLTC we understand a VT-dependent thermal conductivity which is defined 
as follows \v\ 



where 



K (T,V/i;VT) 



< Qx >-- 



l 

V 



T JO 



d<Q x > 
d{V x T) 

T dtQ x (t) 



vivo 



(7) 



(8) 



with Q x (t) being the longitudinal component of the total thermal flux which is defined (in 
a q-space representation) via the total energy conservation law as follows (V is sample's 



volume, and r is a characteristic time of the problem, see below) 

(9) 



Q(t) = lim 

q^O 



with 



(10) 
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Here, s = 2ttc1 2 is properly defined normalization area, and we made a usual substitution 
jf Aij(t) — ► i / cPxA(x.,t) valid in the long-wavelength approximation (q — > 0). 

In view of Eqs.(l)-(10) and assuming 0^- = 7r/2 for the initial phase difference (which 
maximizes the Josephson current of the network), we arrive at the following analytical 
expression for the temperature and chemical pressure dependence of the nonlinear thermal 
conductivity (TC) of a model granular superconductor 



«(T,V/x;VT) 



(11) 



Here, K n = k(T c , V/x; VT) = 2irNd 2 S I c (0)/Vf3c(0) is the normal state value of the NLTC, 
Pc(T) = 2tiIc{T)CqR 2 1 / '$o is the so-called Stewart-McCumber parameter [19j responsible 
for dissipative properties of the array with I C (T) = (2e/h)J(T) being the critical current, 
and 7] = (Efj, — Et)/Eq with Eq = H/2edr being a characteristic field. Notice that, as 
expected, in the limit Et —> (or when VT = in the rhs of Eq.(7)), Eq.(ll) reduces to 
the expression for the linear thermal conductivity k(T, V/x;0). 

Figure [1] shows the dependence of the normalized NLTC k(T, V/x; VT)/k(T, 0; VT) on 
reduced temperature T/Tc for chemical pressure A/x/A/io = 1 (with A/i = h/r) for different 
values of the applied thermal gradient r\ T = St\VT\/E with fi c (0) = 1 and R n = R . Notice 
a markedly different behavior of the nonlinear TC (corresponding to T)t = 0.4, 0.8, and 1.2). 
Unlike its linear counterpart (shown by the upper curve with t]t = 0), it increases with 
increasing the temperature. Even more drastic difference between the linear and nonlinear 
TC can be seen in their chemical pressure dependence. Figure [2] depicts the behavior of 
the NLTC k(T, V/i; VT) as a function of A/x/A/xo for different values of the dimensionless 
parameter 7^ = Et/E and for T = 0.2Tc (the other parameters are the same as before). As 
is clearly seen from this picture, in a sharp contrast with the pressure behavior of the linear 
TC (shown by the curve at the bottom and corresponding to rjT = 0), its nonlinear analog 
evolves with the chemoelectric field quite differently. Namely, NLTC exhibits chemically 
stimulated enhancement with a pronounced maximum at E^ = A/x/2ed ~ Et- To complete 
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our study, let us estimate the order of magnitude of the main model parameters. Starting 
with chemoelectric fields needed to observe the above-predicted nonlinear field effects in 
granular superconductors, we notice that according to Figure [21 the most interesting behavior 
of NLTC takes place for ~ E . Using typical YBCO parameters 1(3], e(0) = 0.01, 
Qo = a|j with ao = 0.2nm, and K = 115GPa, we have fi v = e(0)KQ ~ lmeV for an estimate 
of the chemical potential in HTS crystals, which defines the characteristic time r ~ h/fi v ~ 
5x 10~ n s. Furthermore, taking d ~ lOnm for typical values of the average grain size (created 
by oxygen-rich superconducting regions), we get E = h/2edr ~ 5 x 10 5 V/m and |V/i| = 
Afi/d ~ [i v /d ~ 10 6 eV/m for the estimates of the characteristic field and chemical potential 
gradient (intrinsic chemical pressure), respectively. On the other hand, the maximum of 
NLTC occurs when this field nearly perfectly matches an "intrinsic" thermoelectric field 
Et = jSn|VT| induced by an applied thermal gradient, that is when ~ E ~ Et- 
Using 20j Sq ~ 0.5[iV/K for an estimate of the linear Seebeck coefficient in YBCO, we 
obtain |VT| ~ Eq/Sq — 2 x 10 6 Kjm for the characteristic value of applied thermal gradient 
needed to observe the predicted here effects. Let us estimate now the absolute value of 
the normal state thermal conductivity K n = 2irNd 2 SqIc{0) /V f3c{0) ■ Recall that within the 
present scenario, the scattering of normal electrons is governed by the Stewart-McCumber 
parameter Pc(T) = 27r/c(T)Co-R^/$o ^ ue ^° ^ ne presence of the normal resistance R n 
and mutual capacitance Cq between the adjacent grains. The latter is estimated to be 
Co — laE using d = lOnm for an average "grain" size. Furthermore, the critical current 
Jc(0) can be estimated via the critical temperature Tq as follows, /c(0) — 27rfcsTc/$o 
which gives /c(0) ^ 10/iA (for T c ^ 90K) and leads to /5c(0) ~ 3 for the value of the 
Stewart-McCumber parameter assuming R n ~ R for the normal resistance which, in turn, 
results in q ~ & /R n ~ 10~ 19 C and Ec = q 2 /2Cq ~ O.leU for an estimate of the "grain" 
charge and the Coulomb energy. Finally, assuming V ~ Nd s for the sample's volume, and 
using the above-mentioned expressions for So and /3c(0), we obtain n n ~ 10~ 3 W/mK for 
an estimate of the maximum of the NLTC which is actually much higher than a similar 
estimate obtained for inductance controlled n n in electric- field driven NLTC [171 ]. suggesting 
thus quite a realistic possibility to observe the predicted here non-trivial behavior of the 
thermal conductivity in non-stoichiometric high-Tc superconductors. 
This work was financially supported by the Brazilian agency CAPES. 
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FIG. 1: The dependence of the nonlinear thermal conductivity on reduced temperature T/Tc 
for chemical pressure Afi/AfiQ = 1.0 and for different values of the applied thermal gradient 
rjT = Sq\VT\/Eo (t]t = 0,0.4,0.8, and 1.2, increasing from top to bottom), according to Eq.(ll). 
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FIG. 2: The dependence of the nonlinear thermal conductivity on the chemical pressure A/u/A^o 
for different values of the applied thermal gradient t\t = So\VT\/Eq {ryr = 0,0.4,0.8, and 1.2, 
increasing from bottom to top), according to Eq.(ll). 
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